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We theoretically analyze the temperature behavior of paraxial light in thermal 
equilibrium with a dye-filled optical microcavity. At low temperatures the photon 
gas undergoes Bose-Einstein condensation (BEC), and the photon number in the 
cavity ground state becomes macroscopic with respect to the total photon number. 
Owing to a grandcanonical excitation exchange between the photon gas and the 
dye molecule reservoir, a regime with unusually large fluctuations of the condensate 
number is predicted for this system that is not observed in present atomic physics 
BEC experiments. 



For many problems in statistical physics, one has the freedom of choice to use different 
statistical ensembles for their description, as they often predict the same physical behavior 
in the thermodynamic limit. An interesting counterexample (among others) is Bose-Einstein 
condensation |2| where the grandcanonical description leads to unusually large condensate 



number fluctuations of order of the total particle number 3^6| - in contrast to the predictions in 
the (micro-) canonical case that is typically realized in experiments with ultracold atomic Bose 
gases [7]. This peculiar behavior is known as the grandcanonical fluctuation catastrophe J|3, 6|. 
In recent work, we have observed Bose-Einstein condensation of a two-dimensional photon gas 
in an optical microcavity 

as. 

Here, the transversal motional degrees of freedom of the 
photons are thermally coupled to the cavity environment by multiple absorption-fluorescence 
cycles in a dye medium, with the latter serving both as a heat bath and a particle reservoir. 
Due to particle exchange between the photon gas and the molecular reservoir, grandcanonical 
experimental conditions are expected to be realized in this system. 

In this Letter, we discuss the thermalization mechanism and derive statistical properties 
of the photon condensate, including its photon number distribution, fluctuations and inten- 
sity correlations. The main result is that photonic Bose-Einstein condensates, owing to the 
grandcanonical nature of the light-matter thermalization, can show unusually large particle 
number fluctuations, which are not observed in present atomic Bose-Einstein condensates. 
Our calculations are done in the limit of a non-interacting photon gas. 
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Figure 1: (a) Optical microresonator enclosing a dye medium (left). For a fixed longitudinal mode 
number q, the photon gas can be mapped onto a two-dimensional gas of massive particles (effective 
mass m) confined in a harmonic trap with trapping frequency £1 (right) (b) Scheme of the density of 
states in the resonator. 

The system under investigation, as shown in Fig. la, consists of a microresonator formed 
by two perfectly reflecting spherically curved mirrors enclosing a dye medium. In the cavity, 
optical photons are permanently absorbed and re-emitted by dye molecules. For a sufficiently 
small mirror spacin g, t his process will maintain the longitudinal mode number of the photons 
denoted by q, see [8-10 1 for a detailed description of the experiment. The photon gas effectively 
becomes two-dimensional, as only the two transversal degrees of freedom remain. Throughout 
this Letter, we do not consider experimental imperfections as mirror losses or non-radiative 
decay of dye excitations. Moreover, we do not consider coherent time evolution of the combined 
light-dye system (strong coupling) which would result in polaritonic eigenstates. Under the 
experimental conditions of 8|, |9|, this is inhibited by a rapid, collision-induced decoherence 
process. 

There is a close connection to the statistical physics of the atomic two-dimensional Bose 
gas, which can be seen by expressing the photon energy as a function of the longitudinal (k z ) 
and transversal wave number (k r ), E=hcyjk 2 +k 2 ) with c as the speed of light in the medium. 
The boundary condition in the z-direction is incorporated by an ansatz k z (r)=qn/ D(r), where 



D(r) = D —2(R—\/ R 2 —r 2 ) is the mirror separation at a distance r from the optical axis and 
R is the radius of curvature. In a paraxial approximation, with k r <C k z (0) and r ^ R, this 
yields H 



E 



n (flk 

mc + 



+ 



(i) 



2m 2 

which is formally equivalent to the energy of a two-dimensional harmonic oscillator when 



defining a photon mass m = hk z (0)/c and a trapping frequency Q = c/ ' \J DqR/2. The photon 
gas can thus be mapped onto a (2d) gas of massive bosons in a harmonic trap - a system 
known to undergo BEC at a non-zero temperature 

HQ 



Thermal equilibrium and condensation. - First we discuss the thermalization process 
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in more detail. Thermal equilibrium between photons and dye molecules will be shown 
to arise under two conditions: (i) the Einstein coefficients of the dye medium fulfill the 
Kennard-Stepanov law and (ii) chemical equilibrium between photons, excited and ground 
state molecules is achieved. 

(i) The Kennard-Stepanov law [13I (and references therein) relates the Einstein coeffi- 
cients of stimulated absorption and emission, and can be stated in the form 

■#H =— e kBT > 2 

B 12 {u) w t 

where -612,21 (^) are the Einstein coefficients of absorption and emission, Uq is the zero-phonon 
line of the dye and = f e>0 T>^(e) exp (— e/k^T) de are statistical weights related to the 
rovibronic density of states T>^(e) of ground (\) and excited (7) dye state. This law is often 
fulfilled for dyes in liquid solution and goes back to a thermalization process of the rovibronic 
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dye state due to frequent collisions with solvent molecules 

(ii) The excitation exchange between photon gas and molecules can be seen as a photo- 
chemical reaction of the type 7+ \, ±=t 7. Here 7 stands for a photon, 4- for a ground state 
molecule, and 7 for an excited molecule. The energy of optical photons is well above thermal 
energy at room temperature, i.e. E 3> HbT, and the conversion of ground state molecules 
into excited molecules by thermal fluctuations is thus negligible. In this situation, the total 
number of photons and excited molecules is not adjusted by temperature as for black-body 
radiation, but is a conserved quantity (it depends on the initial state of the photon gas and 
the molecular medium). The chemical potential of the photons /i 7 then can become non-zero. 



In equilibrium, /x 7 is related to the other chemical potentials by /i 7 + //j, = p^ jl6]. Starting 
from this, the photon chemical potential /i 7 can be shown jl5| to be related to the excitation 
level p-j- / p| in the dye medium by 

&r W i Pt P% f \ 
e k B T — __L _J_ e k B T _ (3j 

Pi 

Here are the densities of ground and excited state molecules. The chemical potential p 1 
is spatially homogeneous in equilibrium (as is the temperature), from which follows that the 
excitation level of the medium p\/ p\, is also position-independent. 

The thermalization process is considered as a random walk in the configuration space 
of all allowed light field states. Here, a state K is given by the cavity mode occupation 
numbers K = (riff, nf, , ■ ■ ■)■ The mode occupation numbers are frequently altered by 
photon absorption and emission processes. In first order perturbation theory, the rates (per 
volume) for absorption and emission of one photon in mode i at cavity position r, denoted by 
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i?^' l (r) and R 2 {\r), have the form 

itf/(r) = B^A^Pinf (4) 
</(r) = 5 21 (^K(r)p t (nf + 1) , (5) 

where u l {r) is the spectral energy density of one photon in mode i. We first assume a grand- 
canonical ensemble limit, i.e. we consider the number of dye molecules as sufficiently large, 
that the change of the excitation level p^/pi from photon absorption and emission can be 
neglected. Thus, the species densities p^ are treated as fixed parameters. 

In the theory of random walks, a well known detailed balance criterion exists that deter- 
mines if the rates given by eq. (J3J) and ([5]) lead to equilibrium [l?} Suppose that a state 
K' emerges from state K by the absorption of a photon in mode i, with nf = nf — 1, and 
accordingly that K emerges from K' by an emission process into this mode. The correspond- 
ing (local) rates are Rf 2 {r) and R 21 ' l (r); and with eq. (jlj) and their ratio is given by 
-^i2 ,l ( r )/-^2i '*( r ) = Bni^i) Pi/ B 2 \(uJi) p-|-. Thermal equilibrium will be reached, if this ratio is 
given by the Boltzmann factor of the energy difference between K and K' [18], i.e. 

B 12 (uJi) Pi 



(6) 



B 2 i(uJi) p t 

If one now applies the Kennard-Stepanov relation, eq. fl2]), and assumes chemical equilibrium, 
eq. ([3]), the detailed balance condition eq. ([6]) is indeed verified. Thus, the state of the 
photon gas will thermalize to the temperature of the dye solution T at a chemical potential 
p 7 related to the excitation level of the dye molecules. In particular, the average occupation 
number of mode i, denoted by Hi, can be determined by balancing the average absorption and 
emission rates at a given cavity position. This gives the expected Bose-Einstein distribution 
fii = (exp [{hjJi — p 7 )//cbT] — l) . For atomic Bose gases, one typically omits the rest energy 
mc 2 from the Hamiltonian. We will adopt this convention for the photon gas by removing the 
ground state energy E q oo from all energies including the chemical potentials. With p = p 7 — E q oo, 
the Bose-Einstein distribution becomes 

n{u) = , (7) 

e k B T _ I 

where u = E — E q oo = 0, M7, 2hfl, ... is the reduced photon energy, and we have used the 
degeneracy g(u) = u/hQ + l to combine all modes of same energy, see Fig. lb. For the sake of 
simplicity, we neglect the polarization degeneracy of each mode throughout this Letter, which 
would give an additional degeneracy factor of 2. This case is discussed in the Supplemental 
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15]. For this system, a phase transition is expected at a critical temperature of 
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at which the ground state occupation uq = n(0) becomes a macroscopic fraction of the total 
average photon number N. 

Condensate fluctuations. - For grandcanonical Bose-Einstein condensation, i.e. in the 
presence of an infinitely large particle reservoir, condensate number fluctuations of order of 
the total particle number occur {si . However, no excess fluctuations occur in the (micro- 
canonical ensemble, which is usually realized in atomic BEC experiments [7]. For the photon 
Bose-Einstein condensate, the ground state mode is coupled to the electronic transitions of 
a given number of dye molecules. In this way, the condensate exchanges excitations with a 
reservoir of size M, which is given by the product of dye concentration and ground state mode 
volume. 

We start with the master equation for the probability p n = p n (t) to find n(= no) photons in 
the ground state at time t. The flow of probability between the ground mode and its reservoir 
is governed by [19| 

Pn = Rt-tPn-1 ~ {Rl?+K)Pn + <+lPn+l , (9) 

where i?* 2 = B 12 (M —X +n)n is the rate of absorption and R 2 } = B 2 \{X— n)(n + l) the rate of 
stimulated and spontaneous emission for a configuration with n photons, X—n electronically 
excited molecules and M— X+n ground state molecules. The excitation number X is the sum of 



the ground mode photon number and t 
calculation, X is treated as a constant 



re number of excited molecules in the reservoir. In this 



15] . i.e. it is not expected to perform large fluctuations 



on its own. The rates R} 2 and R 2 } follow from eq. (jl]), © by integrating over volume, 
and setting B 12t21 : = B 12i21 (E qQ0 /h) u q00 (0), where u q00 (0) is the spectral energy density of the 
ground mode on the optical axis. For large times, p n (t) is expected to become stationary, 
p n (oo) = 0, and approach its equilibrium value V n :=p n (oo). In this asymptotic case, the 
solution of the master equation eq. (J9j) is found to be V n = Vq rifc=o •^fV-^ifc+i) wn i cn is 



Vn (m-x)\x\ ihi 

V " (M-X+n)\ {X-n)\ { B, 



(10) 



This photon number distribution can be used to obtain both the average ground state oc- 
cupation and its fluctuations. We consider experimental conditions, where the temperature 
of the system is varied, with the total average photon number N (in all modes) being fixed. 
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Figure 2: (a) Condensate fraction Tlq/N (solid lines) and fraction uq/X (dashed lines) as a function 
of the reduced temperature T/T c for a photon number of N = 10 , six reservoir sizes from Mi = 10 8 , 
M2 = 10 9 up to Mq = 10 13 and a dye-cavity detuning of h/S. / k^,T c = — 4.35. (b) Normalized condensate 
fluctuations (0) = (no (no — 1)) /n§ as a function of T/T c for various M (as in (a)), revealing large 
number fluctuations even far below T c . 
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Figure 3: Phase diagram of the two-dimensional photon gas for fixed average photon number N, as 
a function of the reduced temperature T/T c and the dye-cavity detuning hA/k^Tc. The solid line 
marks the BEC phase transition. The dashed lines (three cases are shown) separate two regimes: 
a condensate with large number fluctuations and a BE-like photon number distribution C(I), and a 
non-fluctuating condensate obeying Poisson statistics C(II). The temperature of the crossover C(I)- 
C(II) depends on the ratio \JM /N , with the reservoir size M as the number of dye molecules in the 
mode volume of the ground state. The insets give a sketch of the corresponding temporal evolution 
of the condensate photon number no(i). 

Correspondingly, upon temperature variations the excitation level p^/pi of the dye medium 
and the total ground mode excitation parameter X have to be readjusted (see Supplemental 



Material 



15]). 



Figure 2a shows the ground mode occupation Uq/N (solid lines) as a function of the reduced 
temperature T/T c for a constant average photon number of N = 10 4 , six reservoir sizes from 
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Mi = 10 8 over M 2 = 10 9 up to M 6 = 10 13 , and a dye-cavity detuning of HA/k B T c := (-EUp- 
huuo) I 'kftT c = —4.35, which adjusts B 2 i/Bi 2 via eq. [2] and is experimentally realizable ja, [9|. 
Upon reaching criticality, the ground mode occupation becomes a macroscopic fraction of 
the total photon number. The occupation level closely follows the analytic solution uq/N — 
1 — (T/T c ) 2 (not shown in Fig. 2a), nearly independent of the reservoir size M (the solid lines 
for different M essentially overlay in Fig. 2a). The dashed lines shows the fraction Uq/X 
of ground mode occupation and total excitation number versus temperature, which reveals 
that the vast majority of excitations is stored as electronic excitations in the medium down to 
relatively low temperatures. Figure 2b gives the normalized ground mode fluctuations in terms 
of the zero-delay autocorrelation function <^ 2 )(0) = (no (no — l))/ng versus temperature. For 
T>T C , one finds the usual case of strong intensity fluctuations with g^(0) =2, accompanied 
by a Bose-Einstein-like photon number distribution. Interestingly, for large reservoir sizes M 
the intensity correlation function remains at this value even at temperatures below T c , i.e. 
when condensation sets in. In this regime, large condensate number fluctuations occur due to 
grandcanonical particle exchange with the molecular reservoir. At even lower temperatures, 
T <C T c , the fluctuations are damped, and one finds <7^(0) ~ 1, accompanied by a Poissonian 
photon statistics. We have no indication that the crossover from Bose-Einstein to Poissonian 
statistics is accompanied by an additional phase transition, i.e. a non-analytic behavior of 
thermodynamic functions. 

In general, the photon number distribution eq. ffTO]) can either be Bose-Einstein-like, 
Poisson-like or an intermediate case. It will be Bose-Einstein-like, if the molecular reser- 
voir is so large that the ratio V n +i/V n = (-B21/-B12) (X — n)/ (M — X + n + l) is approximately 
constant for all relevant photon numbers n. Then V n has its peak value at n = and decays 
exponentially like a geometric series, with V n /Vo ~ [(-B21/-B12) X/(M— X)] n . Thus, the most 
probable event is to find no photons at all. Note that this can be achieved by increasing 
M, while keeping the excitation level Pf/pi — X/(M — X) fixed (which maintains fi, Uq, N). 
For Poissonian statistics, V n has its maximum at a non-zero n. Due to the smooth crossover 
the distinction between these two statistical regimes is not unambiguous. However, a natural 
choice is to consider the point at which 'finding zero photons' ceases to be the most probable 
event. This occurs at Vq = V\ (compare also to the laser threshold described in [20|) and with 
eq. (llOp corresponds to (M + l)/X = 1 + B 2 i/Bi 2 . The temperature T x at which this occurs, 
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for fixed TV and M, can be shown to fulfill the transcendental equation 



6 \ Kl J y ^ + e fc B T x ^x + e vx) 

For zero dye-cavity detuning A = 0, one finds the analytic solution T x a^=o — T c J l — y/M/2N, 
provided that \/ r M/2N < 1, while for general detunings A, eq. (ITT]) has to be solved 
numerically. Figure 3 gives a phase diagram, where solutions for three different cases 
y/~M/2N = \/0.1, 1, x/IO are marked as dashed lines, which separate two regimes of the 
condensate, denoted by C(I) with Bose- Einstein-like photon statistics and C(II) with Pois- 
son statistics, respectively. In terms of second order correlations, the dashed lines corre- 
spond to # (2) (0) ~ 1.571. Note that both T c and T x are conserved in a thermodynamic limit 
N, M, R— > oo that includes N/R = const and yM/N = const. For a fluctuating condensate 
in the C(I) regime, the time dependent intensity correlations g^ir) = (no(t) (no(t+r) — l)) /t%q 
decay like g^(r) = l+exp(— t/t c ), with =no/ B21X as the second order correlation time 
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[15], similar to the photon bunching of thermal emitters. 

To conclude we have studied the thermalization and condensation of photons in equilibrium 
with a dye microcavity. Our analysis predicts a regime with unusually large condensate fluctu- 
ations, not observed in present atomic BEC experiments. Moreover, it adds to a classification 
of photonic BEC in relation to other light sources. In general, Bose-Einstein condensation 
applies to a system in thermodynamic equilibrium; in contrast to lasing which normally oc- 
curs under non-equilibrium conditions. The latter is brought about by a violation of chemical 
equilibrium which occurs if photons are lost (e.g. transmitted) instead of being reabsorbed by 
the medium. Additionally, BEC shows a pronounced dependence on the density of states of a 
system intrinsically requiring a statistical multi-mode treatment, whereas single-mode treat- 
ments are fully sufficient to describe the (single-mode) behavior of a wide class of lasers. If 
one considers a zero-delay autocorrelation function of g^ 2 \0) = 1 to be an essential ingredient 



of lasing, as it is common in laser literature 



2lfl , the intensity fluctuations (photon statistics) 



can give a further distinguishing feature between lasing and photonic BEC, as the latter is not 
necessarily accompanied by a damping of fluctuations. We note that for other lasing threshold 
definitions based on the internal operation of the device, also laser light with non-poissonian 



photon statistics can occur 
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23]. This ambiguity does not occur in the case of photonic 
BEC, where the phase transition is uniquely defined. 

For the future, it will be important to study the statistics of an interacting photon gas. In 



general, interactions can suppress above-poissonian intensity fluctuations [6]. One here has to 
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distinguish between thermooptically induced interactions 8], which are comparatively slow, 
and ultrafast Kerr-lensing. Another line of research could also be dedicated to the (interaction 
induced) superfluidity of the photon condensate. 

We acknowledge financial support of the DFG under contract Wel748/17. 
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